Muonic hydrogen cascade time and lifetime of the short-hved 25" state 
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Metastable 2S muonic-hydrogen atoms undergo coUisional 25'-quenching, with rates which depend 
strongly on whether the fip kinetic energy is above or below the 2S — » 2P energy threshold. Above 
threshold, coUisional 25* 2P excitation followed by fast radiative 2P —r IS deexcitation is allowed. 
The corresponding short-lived fj,p{2S) component was measured at 0.6 hPa H2 room temperature 
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gas pressure, with lifetime rlg""^ = 165_29 
density) and population efs"* = 1.70+° **° 
time, Tcl^ = (37 ± 5) ns, was found. 
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Muonic hydrogen {fi~p) is a simple atomic system, sen- 
sitive to basic features of the electromagnetic and weak 
interactions. Of particular interest is its metastable 25* 
state, long sought after for measuring the /ip(2S')-Lamb 
shift, Cfj^. Vacuum polarization dominates It shifts 
the 25* level by -0.2 eV below the 2P level 0,0. A 
measurement of C^, is in progress at the Paul Scherrer 
Institute (PSI), Switzerland [3|. We report here the data 
analysis of a preliminary stage of this experiment made 
at low II2 gas pressure pu^ =0.6 hPa 

Muons stopped in II2 gas form highly excited fip 
atoms (H. A cascade of both coUisionally- induced and 
radiative deexcitations leads to the IS' ground state or, 
with a probability 623 (few %), to the 25 state. The 
lip{lS) kinetic energy distribution E^f^^ has been mea- 
sured at = 0-06 • • • 16 hPa Q, and cascade simula- 
tions show that E^f^ and E^f^ do not differ significantly 
under our conditions 

The 25* state lifetime is, in absence of collisions, es- 
sentially equal to the muon lifetime t^^ — 2.2 fis. In II2 
gas, coUisional 2S'-quenching occurs, with different pro- 
cesses for kinetic energies E^f^-^ above or below the 25* — > 
2P transition threshold which is (1 -f m^p/mii^)\C^\ ~ 
0.3 eV in the lab frame: 

(i) Most ^p(2S) atoms are formed at energies above 
this threshold [6| where a coUisional 2S 2P Stark tran- 
sition, followed by 2P IS deexcitation with 1.9 keV 
Ka X ray emission (the 2P-lifetime is 8.5 ps), 

fipC^S) + H2 ^ pp{2P) + H2 ^ fipi^S) + H2 + Ka , (1) 



leads to fast 25-dcpletion (coUisional quenching). A life- 
time t|^°''' ^ 100 ns/pH2JhPa'] was predicted for this 
short-lived 2S'-component 0, [1, @1- This is too short to 
have been seen in previous searches for Ka x rays delayed 
with respect to the /ip(2S') formation time 13, In 
this Letter we report on the first measurement of T|g°''* 
and the corresponding population £25™*- 

(ii) Due to clastic collisions, a fraction of the fip{2S) 
atoms decelerates to energies below 0.3 eV where pro- 
cess ^ is energetically forbidden. This fraction is 
the long-lived 2S'-component. A recent experiment [l3l | 
showed that its dominant quenching process is nonra- 
diative deexcitation to the ground state, with lifetime 
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The cascade time T^tP is the mean delay between /ip- 
atom formation and final deexcitation to the ground state 
(when a fip /\-linc x ray, other than from fip(2S) decay, 
is emitted). The T^p. value results from the average of 
the various cascade deexcitation processes, and depends 
on • It was calculated ^ but never measured for fip. 

In our experiment, muons stop in II2 gas containing a 
small admixture of N2 (air), and we measure simultane- 
ously the three time distributions: 

- Ka , i.e., fj.p{n = 2^1) X rays (1.898 keV), 

- i<>;3, i-e., fipin > 3 1) X rays (2.45(2) keVEi), 

- fj,N , i.e., fiN{n = 5 ^ 4) x rays (3.08 keV [1^ 

The fip{n = 3^1) Ki3-\me (2.249 keV) is not weU sep- 
arated from Ka and K^p. The 0.4(1)% air admixture 
in the H2 was useful for calibrating x-ray energies and 
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times. The time distribution is similar to that of 
/ip formation because the fiN cascade time is negligibly 
short 10^^° s) The ^ip cascade time will therefore 
show up as a time delay in the and if>^ distributions 
compared to fiN. The signature for 25 radiative decay 
will be a tail in Ka not present in A>/3. The muon trans- 
fer rates /xp + N — > + p are ~ 10'^ for fip{lS) and 
-lO'^s-i for fip{2S) [i3l, too small to affect our results. 

The experiment was performed at the recently devel- 
oped low-energy ^~ source attached to the ttES beam 
line at PSI [18[. It provides ~ 10^ fj,~ with energies 
of a few keV. The muons were axially injected into a 
1 m long, 20 cm bore solenoid operated at 5 T. con- 
taining the muon entrance detectors and the gas target 
(see [3[). Two detectors, based on nanometer-thick car- 
bon foils, signaled the arrival of slow muons Muons 
were stopped in a 20 cm long target vessel filled with 
0.6 hPa H2 gas (temperature 290 K), and fip atoms were 
formed in a volume of 0.5 x 1.5 x 20 cm*^. Twenty large- 
area avalanche photodiodes (LAAPD), each with sen- 
sitive area 13.5 x 13.5 mm^, were used as x-ray detec- 
tors [2^ . Muon-dccay electrons were detected by a set of 
plastic scintillators and also by the LAAPDs. More than 
5 X 10^ events, each where an x ray was followed by an 
electron, were analyzed. 

Calibration data were used for each LAAPD to deduce 
the energy and time (relative to muon entrance) 
of a measured x ray. Typical resolutions (full width at 
half maximum) were AEx/Ey^ « 25% and At^ ~ 35 ns 
for 2-kcV X rays. Most /jp x rays were found in the time 
interval 300 < tx < 600 ns, corresponding to the widely- 
distributed muon slowing-down times. 

The Ka, Kyfj, and /zN time distributions were deter- 
mined from a fit of the E'x spectra for different tx inter- 
vals. The useful <x range (0.2 • ■ • 6.5 /zs) was divided into 
28 intervals (50, 100, and 500 ns wide). For each interval 
an ii'x spectrum was produced. Three typical spectra are 
shown in Fig. [T] one at times of fj,p formation and deex- 
citation (top), and two at later times. The function fit to 
each spectrum is composed of several x-ray lines and a 
continuous background. Each line is the sum of a Gaus- 
sian peak and a tail towards lower energies (an LAAPD 
characteristic), with energy-dependent weights (3]. 

The lines fitted in Fig. [T] correspond to x rays from fip 
{Ka, Kfj, and K-^p), (main transitions at 1.67, 3.08, 
and 6.65 keV [11]), /iO (2.19, 4.02, and 8.69 keV), and ^iG 
(4.89 keV). The intensities for the Ka, Kyp, 3.08, 4.02, 
4.89, and 6.65 keV lines were free parameters, whereas 
the relative intensities of the other lines as well as the 
positions and widths of all lines were fixed by requiring 
global consistency for all data [3| and considering known 
yields 2l|. The late-time 2-keV peak (see Fig.[Tl bottom) 
is due to ^p X rays from second-muon stops, i.e., muons 
entering the target at random times shortly after a first- 
muon which defined — 0. The ^Cji^s (4.89 kcV) line at 
late times arises from fip atoms drifting to the polypropy- 
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FIG. 1: (color online). X-ray energy spectra (3 of 28) for dif- 
ferent ix-intervals. The fit function is composed of the peaks 
for fip Ka, Kp, K^p, /iN, /iO, /iC(4.89 keV), and a continuous 
background (Bgr). 



lene foils in front of the LAAPDs where muon transfer to 
C atoms occurred. The continuous background is due to 
X rays with energy > 10 keV, e.g., /iC K and L line transi- 
tions, whose deposited charge was not fully amplified by 
the LAAPDs. It was well modeled by a sum of an expo- 
nential and a linear function, with two amplitudes as free 
parameters. The full details of the extensive background 
studies are found in 0|. 

The fitted intensities (with statistical errors) of the 
28 energy spectra are shown in Fig. [2] as a function of tx, 
for the three lines /xN(3.08 keV), fj,p Kyp, and Ka- The 
late-time events are caused by second-muon stops. The 
tx spectra from different LAAPD positions show Q that 
for muons not stopped during their first pass through the 
gas target, a considerable fraction was reflected at the 
gold-plated back side of the vessel and stopped during 
the return pass. Consequently, the fxN spectrum (show- 
ing the muon stop-time distribution) was represented by 
a sum of two functions, each one the convolution of a 
Gaussian with an exponential (Fig. [21 top). The simul- 
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FIG. 2: Fit functions (solid line) and residuals (normalized 
by the errors) for fits to the x-ray time spectra for the lines 
/iN(3.08 keV) (top), /ip — Ky/s (middle), and Ka (bottom). 
Each point results from a fit of the x-ray energy spectrum for 
the corresponding tx-interval (50, 100, and 500 ns wide). The 
dotted lines in the fiN spectrum are the functions for incoming 
(inc), reflected (ref), and second (bg) muons. In the _K>^ 
spectrum the dotted line is the fit function without cascade 
time. In the Ka spectrum the 2S-tail, i.e., the contribution 
of the short-lived ^,p{2S) state, is shown as dotted line. The 
dashed line is the fit function minus the 2S-tail. 



tancous fit of the three time spectra, using the common 
muon stop-time distribution, gives values for the intensi- 
ties A^N, A>/3, and A^. 

The Kyp and Ka spectra are delayed and sUghtly 
broadened with respect to /xN due to the np cascade time. 
It is not possible to extract the precise shape of the cas- 
cade time distribution from the data, but calculations 0| 
indicate that the K-^p cascade has an approximately ex- 
ponential time distribution, whereas the Ka cascade has 
the same asymptotic behaviour but includes a "build-up" 
character at earlier times. The /iN fit function was there- 
fore convoluted with an exponential (parameter Teas) to 
obtain the Kyp function and further convoluted with a 
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FIG. 3: Relative population £|g°'* (normalized to Aa) of the 
short-lived ^p{2S) component versus its lifetime f||°'' (with- 
out systematic corrections). The dot represents the best fit. 



Gaussian (parameter cr^as) foi' tl^c Ka function. In addi- 
tion, free time offsets Ar>/3 and Ar„ (for K^p and A'q, 
with respect to /iN) were introduced. The contribution 
of the short-lived fip{2S) state was considered by adding 
to the Ka function a convolution of the K^p distribution 
with an exponential (parameter "rl^"'*), and a relative 
population e^^™* (normalized for the fit to Aa). 

A simultaneous fit of the three time spectra was per- 
formed, resulting in ~ 57.5 for 67 degrees of free- 
dom. The fit functions and residuals are shown in 
Fig. [21 The scatter of the residuals confirms that no 
relevant systematic deviations exist between the data 
and the model. The resulting cascade time slope is 
Teas = (26 ± 2) ns. The time shift AT>/3 = (0 ± 5) ns 
is consistent with zero, as expected for the K^p cascade, 
whereas ATq = (13 ± 5) ns and = (15 ± 2) ns ap- 
proximate a Ka cascade time distribution which deviates 
from an exponential in the first ~ 20 ns. The deduced 
mean cascade times Tjas = '^cas + Ar>/3 = (26 ± 5) ns 
and T;^^^ = Teas + ATa = (39 ± 5) ns, weighted by the 
corresponding if -line yields, result in a mean fip cas- 
cade time r^ls = (37 ± 5) ns. (The relative if-yields at 
= 0.6 hPa were deduced from [li| as = 0.821(12), 
Yp = 0.061(9), and Y^p = 0.118(11), and the A> cascade 
time was assumed to be equal to T^.) T^'^p depends only 
weakly on the fit function details because it essentially 
reproduces the center-of-gravity shifts of the Ka and Kyp 
distributions relative to fiN. 



The fit results for the 2S'-tail are 



133+22 ns 

and efs°'^ = 1.44l° j^%. As shown in Fig. [3 the absence 
of the 2S'-tail, corresponding to e^g"^' = 0, is excluded by 
6(T. As a test, the Kyp spectrum was also fit with a "25"'- 
tail. Its amplitude (normalized to Ayp) is compatible 
with zero (within 0.7(t) for any t||°''*, as expected. We 
conclude that our fit function correctly reproduces the 
muon stop, cascade, and second-muon time distributions. 
The difference between e|5°''* and e>;3 is e^'j"'' - e^p = 
(1.36 ± 0.28)% at tI^°'^ = 133 ns. A zero value of this 
difference is excluded by > 4.4 ct for any r|g°''* , confirming 
that the tail in the Ka spectrum can only come from 25- 
quenching. 

The measured fl^""^' and e^'g"'' values have to be cor- 
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rected for several effects. The necessary corrections were 
deduced from a Monte Carlo simulation of the experi- 
ment, based on the known distribution of E'^f^ @ and 
on calculated cross sections 0] for process ([T]) and elastic 
collisions. It was found that (i) some fj,p(2S) atoms reach 
the target walls before being quenched; (ii) the solid an- 
gle for Ka detection varies with time due to the ^p{2S) 
motion; (m) because of collisions, the E^^^ distribution 
depends on time, and hence so do the mean cross sections 
and fip velocities. Those effects modify the 2S tail shape 
at the late times where the experiment is most sensitive. 
The resulting correction factors arc 1.24 ± 0.07 for f||°''* 



and 1.34 ± 0.09 for e: 
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In addition, 



ef^'^ has to be 



multiphed by (1 — T|g /r^) to account for muon decay 
and hy Ya / to normalize to all fip atoms. The fi- 
nal result for the radiative lifetime and population of the 
short-lived fip{2S) component at 0.6 hPa (temperature 
290 K) is t||°'' = I65tll ns and e^^"* = l-70tj] %. 

The fip cascade time and the 2S'-tail have been disen- 
tangled from the stop-time distribution for the first time. 
This has been made possible by the high statistics, the 
low gas pressure, and the good stop-time resolution. The 
measured jip cascade time T/fP^ = (37 ± 5) ns is less than 
half the « 90 ns value predicted for 0.6 hPa by cascade 
calculations 0]. The measured value may be slightly af- 
fected by a possible difference in the atomic capture 
times of N2 and II2 predicted by some models 
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muon energy from which capture occurs is expected to in- 
crease with Z [23! , so /^N atoms can form earlier than ^p 
atoms during the muon stopping, an effect which would 
result in an even lower measured T^^^ value. A result 
which does not depend on such effects is the measured 
difference — = (13 ± 4) ns, also significantly 
smaller than the calculated value of w 25 ns. The calcu- 
lated cascade times may be too long since Coulomb de- 
excitations, which dominate the cascade at high n-levels 
even at low p^.^ , may be accompanied by simultaneous 
Auger transitions, an effect not yet considered. 

We have compared the measured radiative 2S'-lifetime 
with the results of calculations. Since the existing cal- 
culated cross sections neglect molecular effects we al- 
lowed for the extreme cases where the H2 molecule is 
taken as two separate H atoms or as a single atom. 
Analysing the simulated data in the time region where 
the experiment is sensitive, we obtained a value of 
(178 ± 30) ns which agrees well with the experimental re- 



sult 
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= 1651 29 ns. This confirms the validity of the 



cross sections calculated for iip{2S) -I- H — > fj,p{2P) + H 
Stark transitions [1, Q. The quenching rate for the 
short-lived 2S'-component, when normalized to liquid- 
hydrogen atom density (LHD, 4.25 x 10^^ atoms/cm^). 



is A5^°"=^LHD) 
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20 times faster 



than for the long- lived one [13l |. 

The sum of the measured relative populations e^ls"^* ~ 
l-70to:56% and e^°s^ w 1% (extrapolated to 0.6 hPa, 



from [13|) is (2.7 ± 0.8)%, in agreement with the total 
2S'-population £25 = (2.49 ± 0.17)% at 0.6 hPa deduced 
directly from the measured /ip i^T- line yields 14 1 . 

We conclude that there is now good understanding of 
both the short- and long-lived ^ip{2S) dynamics. 
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